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J/ip production at HERA 
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We consider the J /^l> photo- production data collected at HERA in the light of next-to-leading order 
predictions for the color-singlet yield and polarization. We find that, while the shapes of inclusive 
distributions in the transverse momentum and inelasticity are well reproduced, the experimental 
rates are larger than those given by the color-singlet contribution alone. Furthermore, the next-to- 
leading order calculation predicts the J/i/^'s to be mostly longitudinally polarized at high transverse 
momentum in contrast with the trend of the preliminary data from the ZEUS collaboration. 
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Non-relativistic QCD (NRQCD) provides a rig- 
orous framework to consistently separate high and 
low energy effects in quarkonium production and 
decay The factorization of the different en- 

ergy scales in the effective non-relativistic theory 
is achieved through an expansion in v, the rela- 
tive velocity of the heavy quarks in the quarkonium 
state. If u is a sufficiently small parameter, the non- 
pcrturbativc phenomena can be encoded into a few 
process-independent long-distance matrix elements 
(LDMEs) which can be extracted from experimen- 
tal data. 

Whether such an expansion is well-behaved for 
charmonium, and in particular for the J/ip^ is still 
subject to debate. So far no complete experimen- 
tal evidence for the universality of the LDMEs has 
been reached. In particular, the role of the color- 
octet transitions is not well assessed. Although these 
contributions seem to be required to account for the 
observed Pt spectrum in hadro-production [2, 
they fail to explain the polarization measurement at 
large transverse momentum Q and they do not seem 
to be needed to describe fixed-target experiments Q 
and photo-production [6, 7] data. 

Estimating the impact of the transitions at work 
for J /ip production relies on the accuracy with which 
the corresponding short-distance coefficients can be 
computed. Cross sections at leading order in as 
are normally affected by very large uncertainties 
and cannot give a reliable estimate of the yield. 
In these cases, gathering information on the un- 
derlying production mechanism(s) from data can 
be problematic. The recent computation of next- 
to-leading order (NLO) corrections in J/ij} electro- 
production [l^, [ll| has reduced the previous 
discrepancies between the NRQCD predictions and 
the data. In hadro-production, higher-order as cor- 



rections to the color-singlet production have been 
shown to significantly increase the yield at large 
Pt i3j 3 14 1, giving a better description of both 
the shape and the normalization of the Pt spectrum. 
Remarkably, the color-singlet prediction for the po- 
larization as a function of the Pt also dramatically 
changes once corrections in as are taken into ac- 
count, leading to a better agreement with the data. 

In view of the recent theoretical progress related 
to J/V' hadro-production, it is worth reconsidering 
the photo-production data obtained at HERA. Dif- 
ferential cross sections and polarization observables 
for J / %jj photo-production have been measured both 
by the ZEUS and HI collaborations. Previous com- 
parisons of these measurements with the NRQCD 
predictions [1, 0] suggest that the color-singlet yield 
alone reproduces the experimental differential cross 
sections. However, theoretical uncertainties are too 
large to draw any definite conclusion. 

For the purpose of identifying the transitions that 
dominate J /tp production, it is often useful to an- 
alyze the polarization. Contrary to the differential 
cross sections, one can expect the polarization pre- 
dictions to be less affected by the uncertainties as- 
sociated with the theoretical inputs. It has been ar- 
gued [l^ that the J/ip polarization could be a good 
discriminator between color-singlet and color-octet 
transitions in 7p collisions. However, only leading 
order predictions have been considered so far due to 
the lack of NLO calculations for the polarization for 
either the singlet [3l or the octets Since QCD 
corrections can be very large in some regions of phase 
space, especially for the color-singlet at large trans- 
verse momentum, comparisons with leading order 
predictions are of limited interest. 

As new measurements of the J /ip polarization at 
HERA are still being performed by the ZEUS and 
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HI collaborations, we limit ourselves to a compar- 
ison with the preliminary data [l^. For the color- 
singlet transition, polarization observables in photo- 
production are given here at NLO accuracy in as 
for the first time. 

The relevant kinematic variables for J/iJj photo- 
production are the energy of the photon-proton sys- 
tem (W) and the fraction of the photon energy car- 
ried by the J/ip in the proton rest frame (z), 
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as well as the transverse momentum of the J/tp {Pt)- 
Near the end-point region z ~ 1, Pt — GeV, 
J/V' production is enhanced by diffractive contri- 
butions characterized by the exchange of colorless 
states. These contributions cannot be correctly ac- 
counted for by the NRQCD factorization approach, 
and hence they must be excluded from experimen- 
tal measurements in order to make a meaningful 
comparison. As diffractive production leads to a 
steeper slope for the Pt spectrum of the the 
cut Pt > 1 GeV effectively eliminates most of the 
diffractive events. Whether this cut is enough to se- 
lect a clean sample of inelastic events is, however, 
still the subject of experimental investigation [Tot . 
In the following, we adopt the cuts used by the ZEUS 
collaboration. 

Inelastic contributions to J/^p production include 
the feed-down from excited charmonium states, 
which has not been taken into account in this pa- 
per. Feed-down from ip{2S) is seen to lead to a 15% 
increase in the J/tp cross section 0] , whereas contri- 
butions from the decay of Xc states are expected to 
be smaller (« 1%). The J/ip can also originate from 
a resolved photon. This is not a technical problem 
since it requires a calculation of hadro-production 
that is available [l^. However, these contributions 
arc only significant in the low z region, as we have 
checked by direct computation. Since there is no 
data below z ~ 0.3, we consider only the direct con- 
tributions in this paper. 

The ZEUS collaboration has published a measure- 
ment of the differential cross sections for J/ ip photo- 
production produced from electron-proton collisions, 
with Ep = 820 GeV and Ee = 27.5 GeV 0. In 
their analysis, they consider the range 50 GeV < 
W < 180 GeV, and convert electro-production cross 
sections into photo-production cross sections by di- 
viding by the photon flux integrated in this range. 
This is estimated within the Weizsacker- Williams 
approximation: using Q'""^' — 1 GcV, they obtain 
-Fphoton = 0.0987 for the integrated photon flux. For 
a consistent comparison with this measurement, our 
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TABLE I: Scale dependence of the total cross section 
a{PT > 1 GeV) (expressed in nb). W is set to 100 GeV. 
Other input parameters are explained in the text. 



7P cross sections are convoluted with the photon flux 
using the Weizsacker- Williams approximation, and 
then divided by the integrated flux quoted above. 

The cross section (7(7^ J/ip + X) in the color- 
singlct model [2^ has been computed in Ref. [l3| 
at NLO accuracy in as. We reproduced this com- 
putation by using_^thc method and the results pre- 
sented in Ref. ,12] for the hadro-production case, 
which also allows tracking of the J/ip polarization 
through the angular correlations of the J/ip decay 
products, J/ip . We have compared the re- 

sults for several inclusive observables, finding always 
very good agreement. As a further check and ap- 
plication of our results we also calculated the NLO 
corrections to the decay width of a ^5^' into 7 + 
hadrons. Our result is in excellent agreement with 
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that of Ref. 

In our photo-production analysis, we fix 
{Oji^{^Si[l])) = 1.16 GeV^ and we use the 
CTEQ6M pdf set Given the small range in pT 
available and for the sake of simplicity we prefer 
to use fixed scales. As the central value we choose 
/io = 4mc, i.e., where the sensitivity of the cross sec- 
tion a{PT > 1 GcV) to scale variations is minimal 
(sec Table |l|. Since a large variation is observed 
when we vary the scales in opposite directions, we 
impose the additional condition 0.5 < ^ < 2 in 
our prediction for the differential cross sections. 
We also vary the charm quark mass in the range 
1.4 — 1.6 GeV. Scale and mass uncertainties are 
combined in quadrature. The resulting distributions 
are displayed in Fig. [1] 

For comparison, we also plot the color-singlet pre- 
diction at leading order in as, for which we use the 
CTEQ6L1 pdf set. As can be seen from Fig. [TJ the 
as corrections increase the differential cross section 
in the high Pt region, where the yield is dominated 
by the new channels that open up at order a^. Nev- 
ertheless, the color-singlet yield at NLO clearly un- 
dershoots the ZEUS data. The plots in Fig. [T] dif- 
fer from the comparison presented in Ref. [7[ , where 
rather extreme values for the rcnormalization scale 
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FIG. 1: Differential cross sections in and z, compared 
with the ZEUS data 01- 



were used that have the effect of artificially increas- 
ing the normalization. We do not display here the 
comparison with the and z distributions mea- 
sured by the HI collaboration Q since it shares the 
same features. 

We now turn to the polarization. Experimentally, 
the polarization of the J/i/^'s can be traced back 
by analyzing the angular distribution of the leptons 
originating from the decay of the J/tp. It is con- 
venient to decompose this angular distribution in 
terms of the polar and azimuthal angles 9 and (j) 
in the J/^p rest frame: 



da 
dftdy 



oc 1 + A(y) cos^ 6* + sin 26* cos ( 



sin^ 9 cos 2(f> 



(2) 



where y stands for a certain (set of) variable(s) (ei- 
ther Pt or z in the following). If the polar axis 
coincides with the spin quantization axis, the pa- 
rameters A, /i, v can be related to the spin density 
matrix elements: 
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_ (3) 

The spin information that we extract in this way 
depends on the choice of the quantization axis. Here, 



we decide to work in the target frame {z = ^^^) 
as in the recent analysis performed by the ZEUS 
collaboration. 

As mentioned above, in our NLO computation 
the J/ip decays into leptons. In order to obtain 
A (resp. v) we have integrated the cross sections 
over (j) (resp. 9) and extracted the polarization pa- 
rameters by fitting the resulting distributions. In 
so doing, we have used the same kinematic condi- 
tions as those considered by the ZEUS collabora- 
tion: Ep = 920 GeV, = 27.5 GeV, Pt > 1 GeV, 
z > 0.4 and 50 GeV < < 180 GeV. Notice that 
the cut z < 0.9 previously considered for the dif- 
ferential cross sections has been relaxed here. The 
measurement of the polarization versus Pt is there- 
fore subject to a larger contamination from diflrac- 
tive contributions. 

The NLO predictions of the polarization parame- 
ters associated with color-singlet production arc dis- 
played in Fig. [51 together with the LO predictions 
and the preliminary ZEUS measurements. The band 
comes from the uncertainties associated with the 
choice of the scales - varied in the range defined 
by 0.5/io < p-fjP'r < 2/io and 0.5 < ^ < 2 - which 
is much larger than the mass uncertainty. For some 
specific values of the scales (namely /i^ = 0.5), the A 
and parameters appear to be unphysical in some 
bins. This is due to the fact that in these cases our 
calculation leads to a negative value for the diagonal 
components of the spin density matrix at Pt < 1 
GeV, and hence cannot be trusted in this region. 
Therefore, in the error bands in Fig. [2l we have dis- 
regarded scale choices leading to unphysical predic- 
tions. 

QCD corrections to the color-singlet production 
have a strong impact on the polarization prediction. 
The most spectacular effect comes from the behavior 
of the A parameter at large transverse momentum, 
for which the prediction is rather stable under varia- 
tion of the scales. At leading order in as, the color- 
singlet transition gives a transverse J/-0 at large Pt- 
Once QCD corrections are included, the A parame- 
ter decreases rapidly and has a large negative value 
above Pt = 4—5 GeV. This situation is similar to 
that in hadro-production where the J/^j produced 
via a color-singlet transition is longitudinal at large 

Such a correction 



transverse momentum [Ij 
for the A parameter at moderate and high Pt, as 
well as the decrease at z « 0.8, is not supported by 
the preliminary data from the ZEUS collaboration 
and therefore suggests the presence of other mecha- 
nisms for J/iJj production. In the low z region, the 
scale uncertainty is too large to draw any conclu- 
sion. QCD corrections to color-singlet production 
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FIG. 2: Polarization parameters for color-singlet produc- 
tion at leading order and next-to-leading order in as, 
compared with the ZEUS preliminary measurement. 



also affect the value of the v parameter, which is 
closer to the experimental data in comparison with 
the prediction at leading order. 

In this letter, we have studied J/ip photo- 
production via a color-singlct transition at HERA 
and presented for the first time a comparison with 
the predicted polarization at NLO accuracy. Tak- 
ing into account the as corrections, color-singlet 
production alone does not describe all features of 
the data collected at HERA. With a natural choice 
for the renormalization scale, the predicted rate is 
smaller than data, even though the shapes of the dif- 
ferential distributions are well described. Moreover, 
the preliminary measurement of the J/ip polariza- 
tion by the ZEUS collaboration as a function of the 
Pt shows a very different trend with respect to the 
theoretical predictions. 

New channels appearing at NNLO in ag (includ- 
ing fragmentation processes) may increase the dif- 
ferential cross section, as has been pointed out in 
the hadroproduction case [l^l- Although the kine- 
matics differ in photoproduction, such contributions 
could give an enhancement at large transverse mo- 



mentum {Pt > 5 GeV). However, no reliable esti- 
mate of NNLO contributions in the Pt region acces- 
sible at HERA is presently available. 

The current discrepancies could possibly be solved 
by invoking color-octet transitions, i.e. contribu- 
tions from the intermediate states ^s]^^ and '^-Pf'. 
Unfortunately, any phcnomenological analysis of the 
impact of these contributions on differential cross 
sections and polarization observables is limited by 
the omission of higher-order corrections that are cur- 
rently unknown. A complete a| computation, par- 
ticularly for the prediction of the polarization of the 
J/ip produced via a P— wave color-octet state, would 
be welcome in order to shed further light on the 
mechanisms at work in photo-production. 
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